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We investigated the spin-dependent transport properties of Fe/MgAl2O4/Fe(001) magnetic 
tunneling junctions (MTJs) on the basis of first-principles calculations of the electronic struc- 
tures and the ballistic conductance. The calculated tunneling magnetoresistance (TMR) ra- 
tio of a Fe/MgAl2O4/Fe(001) MTJ was about 160%, which was much smaller than that of a 
Fe/MgO/Fe(001) MTJ (1600%) for the same barrier thickness. However, there was an evanes- 
cent state with Ai symmetry in the energy gap around the Fermi level of normal spinel MgA^CU, 
indicating the possibility of a large TMR in Fe/MgAl 2 O 4 /Fe(001) MTJs. The small TMR ratio 
of the Fe/MgAl2O4/Fe(001) MTJ was due to new conductive channels in the minority spin states 
resulting from a band-folding effect in the two-dimensional (2-D) Brillouin zone of the in-plane wave 
vector (/c||) of the Fe electrode. Since the in-plane cell size of MgAbCU is twice that of the primitive 
in-plane cell size of bcc Fe, the bands in the boundary edges are folded, and minority-spin states 
coupled with the Ai evanescent state in the MgA^CU barrier appear at k\\=0, which reduces the 
TMR ratio of the MTJs significantly. 
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I. INTRODUCTION 



Recent advances in tunneling magnetoresistance 
(TMR) in magnetic tunnel junctions (MTJs) with a 
single-crystal MgO barrier and bcc Fe ferromagnetic 
electrodes make possible the fabrication of ultrahigh- 
speed and high-density magnetic random access mem- 
ory(MRAM) devices^. In Fe/MgO/Fe(001) MTJs, the 
crystal momentum parallel to the layer is conserved be- 
cause of the 2-D periodicity of the system, and the tun- 
neling conductance strongly depends on the symmetry 
of the propagating states in the bcc Fe electrode because 
of the complex band structures of MgO, which causes 
a slow decay of the evanescent Ai state in the barrier 
layer. Furthermore, electrons with an in-plane wave vec- 
tor fey = (0,0) (normal incidence with respect to the 
plane) dominate the tunneling, and bcc Fe has a 100% 
spin polarization in the Ai state at the Fermi level. Thus, 
Fe/MgO/Fe(001) MTJs act as spin filters for the current, 
yielding much larger TMR ratios of over 1000% in their 
ground state 4,5 . 

Achieving a high tunneling magnetoresistive (TMR) 
ratio in MgO-based MTJs requires epitaxial growth of the 
MgO layer with a correct crystalline orientation on ferro- 
magnetic electrodes (bcc Fe and Co^Fei-^). However, it 
has been difficult to grow Fe/MgO/Fe(001) MTJs epitax- 
ially owing to the relatively large lattice mismatch (5%) 
between rock-salt-type MgO and bcc Fe. Recently, an off- 
stoichiometric, normal spinel MgAl 2 04 barrier was grown 
epitaxially on a single-crystal Co2FeAlo.5Sio.5 and bcc Fe 
to explore new materials for the barrier layer of MT Js^i 
Since normal spinel MgAl 2 04 has a lattice constant a of 
8. 16 A, the lattice mismatch with bcc- type ferromagnetic 
metals such as bcc Fe and Co-based full Heusler alloys is 



very small, i.e., less than 1% for a 45° (001) in-plane rota- 
tion, indicating that MgAl 2 04 has the potential to over- 
come the problem associated with MgO barriers. In fact, 
a large room-temperature TMR ratio of over 100% was 
reported in Co 2 FeAlo.5Sio.5/MgAl 2 04/CoFe(001) and 
Fe/MgAl 2 O 4 /Fe(001) MTJs. Furthermore, a relatively 
large bias voltage for one-half the zero-bias TMR ratio 
was obtained at room temperature (RT), which is about 
twice that reported for MgO-based MTJs. 



The band structures of bulk MgAl 2 04 were in- 
vestigated by first-principles density functional 
calculations^, and the results were compared 
with experimental data from vacuum ultraviolet 
measurements^. Normal spinel MgAl 2 04 has an indi- 
rect band gap of about 7.8eV (6.5 eV in calculations), 
and the Mg-0 bond length (1.919A) is shorter than that 
in MgO (2. 102 A), implying a stronger bond, and hence 
a larger band gap. However, the coherent tunneling 
properties of MTJs are still unknown for MgAl 2 04. In 
particular, it is necessary to elucidate the symmetry- 
dependent tunneling through the evanescent states of 
MgAl 2 04 and the physical cause of the relatively large 
TMR ratio obtained in recent experiments involving 
Fe/MgAl 2 O 4 /Fe(001) MTJsi. 



This work aims to determine the coherent tunnel- 
ing properties of MgAl 2 04 and the cause of the large 
TMR ratios in Fe/MgAl 2 O 4 /Fe(001) MTJs. To this 
end, we investigated the spin-dependent transport prop- 
erties of Fe/MgAl 2 O 4 /Fe(001) MTJs on the basis of first- 
principles density functional calculations of the electronic 
structures and the ballistic conductance. 



II. COMPUTATIONAL DETAILS 

We prepared a super cell of a Fe/MgAl 2 O4/Fe(001) 
MTJ containing 11 atomic layers of bcc Fe and 9 atomic 
layers of MgA^C^. The in-plane lattice parameter of the 
super cell was fixed at 5. 733 A, which corresponds to twice 
the lattice constant of bcc Fe (2. 865 A). Since the lattice 
constant of spinel- type MgA^C^O is a=8.16A, the lat- 
tice mismatch between bcc Fe and MgAl204 for a 45° 
(001) in-plane rotation is less than 1.0%. We performed 
first-principles calculations of the super cell using density 
functional theory within the generalized-gradient approx- 
imation for exchange-correlation energy 11 . In order to 
facilitate structure optimization, which is important for 
determining the interface structure, we adopted plane- 
wave basis sets along with the ultrasoft pseudopotential 
method by using the quantum code ESPRESSO 12 . The 
number of k points was taken to be 5x5x1 for all cases, 
and Methfessel-Paxton smearing with a broadening pa- 
rameter of 0.01 Ry was used. The cutoff energies for the 
wavefunction and charge density were set to 30Ry and 
300Ry, respectively. These values are large enough to 
deal with all the elements considered here within the ul- 
trasoft pseudopotential method. 

For the conductance calculations, we considered an 
open quantum system consisting of a scattering region 
having a MgA^CU barrier and junctions with bcc Fe 
attached to left and right semi-infinite electrodes corre- 
sponding to bulk bcc Fe. The conductance was obtained 
by solving the scattering equation with infinite boundary 
conditions in which the wavefunction of the scattering re- 
gion and its derivative were connected to the Bloch states 
of each electrode^. Since our system is repeated peri- 
odically in the xy plane and propagating states can be 
assigned by an in-plane wave vector fen = (k x ,k y ) index, 
different k\\ do not mix and can be treated separately. 
Furthermore, we neglected the spin-orbit interaction and 
noncollinear spin configuration. Thus, we solved the scat- 
tering equations for some fixed fen and spin index on the 
basis of Choi and Ihm's approac h 13 ! 14 . 



III. RESULTS AND DISCUSSION 

First, we investigated the stable structure of 
Fe/MgAl 2 O 4 (001) interfaces. Normal spinel MgAl 2 4 
has two different types of cation sites: tetrahedral Mg 
sites (A-sites) and octahedral Al sites (B-sites). This re- 
sults in two types of termination for Fe/MgAl2 04(001) 
junctions, namely, the A-site (Mg) termination and IB- 
site (Al-O) termination. On the basis of formation 
energy calculations for optimized surface structures of 
MgAl2 04(001), we found that B-site termination is ther- 
modynamically more stable than A-site termination. 
This was also reported for other normal spinel com- 
pounds, such as Fe3O4(001) 15 . Thus, we considered only 
B-site termination for the Fe/MgAl2 04(001) interface. In 
B-site termination, there are four possible configurations, 



(a) Al-top (b) AI-0 hollow (c) AI-AI hollow (d) O-top 






FIG. 1: (Color online) Schematic cross-sectional view of a IB- 
site terminated Fe/MgAb 04(001) interface with (a) Al-top, 
(b) Al-0 hollow, (c) Al-Al hollow and (d) O-top configura- 
tions. 



with the Fe atoms positioned on top of the Al atoms (Al- 
top), Al-Al hollow, Al-0 hollow and the O atoms (O-top) 
of MgAl204. In order to determine the stable interfacial 
configuration at Fe/MgAl2 04(001) junctions, we mini- 
mized the total energy by relaxing all atomic positions ex- 
cept for those in the electrode region by changing the lon- 
gitudinal size of the super cell. Figure 1 shows schematics 
of the Fe/MgAl 2 O 4 /MgO(001) interfacial structure with 
B-site termination. We found that the O-top configura- 
tion is the most stable because of the hybridization be- 
tween the 3d 3z 2_ r 2 orbital of Fe and the p z orbital of O. 
Therefore, in the present study on spin-dependent con- 
ductance and the TMR effect, we adopted the O-top con- 
figuration for the B-site terminated Fe/MgAl2 04(001) in- 
terface. 

Figure 2 shows the local density of states (LDOS) of 
Fe in MgAl 2 O 4 /Fe(001), MgO/Fe(001) and the bulk re- 
gion. The interfacial LDOS of Fe at the interface with 
MgAl204 are modified from those in the bulk region, ow- 
ing to the appearance of the nonbonding Fe-3d x 2_ y 2 and 
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FIG. 2: (Color online) Local density of states (LDOS) of in- 
terfacial Fe at MgAl 2 O 4 /Fe(001) and MgO/Fe(001) interfaces 
together with bulk Fe as a function of energy relative to the 
Fermi energy. 



TABLE I: Calculated majority-spin and minority-spin conductance for parallel and anti-parallel magnetization and TMR ratios 
for Fe/MgAl 2 O 4 /Fe(001) and Fe/MgO/Fe(001) MTJs with the same barrier thickness (~2nm). 



Fe/MgAl 2 O 4 /Fe(001) MTJ 
Fe/MgO/Fe(001) MTJ 



Conductance [S//xm 2 



Parallel magnetization 

Majority spin Minority spin 

0.507 0.0245 

0.113 0.00300 



Anti-parallel magnetization 

Majority spin Minority spin 

0.0994 0.105 

0.00346 0.00324 



TMR ratio 



160 
1630 



Fe-3d xy states around the Fermi level. The spin moment 
of interfacial Fe is 3.08/iB, which is larger than that in the 
bulk region (2.49/iB)- Similar results have been obtained 
for MgO/Fe(001) interfaces^. 

We investigated the tunneling conductance of the 
Fe/MgAl 2 O 4 (~lnm)/Fe(001) MTJ. Figure 3 plots 
the in-plane wave vector (fen) dependence of the 
tunneling conductance at the Fermi level for the 
Fe/MgAl 2 O 4 /Fe(001) MTJ in parallel and anti-parallel 
magnetization configurations. We confirmed in Fig. 3(a) 
that the majority-spin conductance for the parallel mag- 
netization has a broad peak around the center of the 
2-D Brillouin zone. This is the typical behavior of 
the coherent tunneling conductance of Ai electrons at 
fen = (0,0). Furthermore, we found in Fig. 3(b) that the 
minority-spin conductance for the parallel magnetization 
also shows a broad peak at fey =(0,0). This result is 
very different from that for the Fe/MgO/Fe(001) MTJi£, 
where the fey -dependence of the minority-spin conduc- 
tance for the parallel magnetization shows hotspot-like 
peaked structures in the 2-D Brillouin zone with no 
peak at fey =(0,0) because of the absence of the Ai band 
around the Fermi level. Since the tunneling conductance 
in the anti-parallel magnetization shows a combination 
of the features observed in the majority and minority 
spin channels, we obtained a broad peak at fey = (0,0) 
in the fey dependence of the conductance for the anti- 
parallel magnetization as shown in Fig. 3(c). In Table 1, 
we showed the majority- and minority-spin conductance 
and the TMR ratios for Fe/MgAl 2 O 4 /Fe(001) MTJs and 
Fe/MgO/Fe(001) MTJs. The calculated minority-spin 
conductance for the Fe/MgAl 2 O 4 /Fe(001) MTJ for par- 
allel magnetization was one order of magnitude higher 
than that for the Fe/MgO/Fe(001) MTJ. This led to 
TMR ratio of 160% for the Fe/MgAl 2 O 4 /Fe(001) MTJ, 
which was much lower than that for the Fe/MgO/Fe(001) 
MTJ (1600%) with the same barrier thickness. 

Figure 4 shows the imaginary part of the out-of-plane 
wave vector k z (k) at the Fermi levels of bulk MgAl 2 4 
and MgO as a function of the in-plane wave vector 
k x . The positions of the Fermi level in the band gap 
were determined from the LDOS of oxygen atoms in the 
center of the barrier layer of Fe/MgAl 2 O 4 /Fe(001) and 
Fe/MgO/Fe(001) MTJs. Since propagating states in the 
metal electrode couple to evanescent states (k) in the bar- 
rier layer and decay as ~ e~ KZ : , the evanescent state with 
the smallest decay parameter, ft m i n , makes the largest 



contribution to the tunneling conductance in MTJs. We 
found that at (fe x ,fe y )=(0,0), MgAl 2 4 has the smallest 
k with Ai symmetry in the gap connecting the top of 
the Ai valence band to the bottom of the Ai conduc- 
tion band, which is similar to the case for rock-salt-type 
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(c) Anti-parallel magnetization 
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FIG. 3: (Color online) In-plane wave vector k\\—(k x ,k y ) de- 
pendence of (a) majority- and (b) minority-spin conductance 
at the Fermi level for Fe/MgAl 2 O 4 /Fe(001) MTJs with par- 
allel magnetization and (c) anti-parallel magnetization. 




(a) majority spin 



(b) minority spin 



0.0 0.1 0.2 0.3 0.4 0.5 

In-plane wave vector k x [2n/a] 

FIG. 4: (Color online) Imaginary part of k z as a function 
of real k x at the Fermi levels for bulk MgO and MgAbO^ 
Note that the Fermi levels for bulk MgO and MgAbCU are 
determined from the LDOS of oxygen atoms in the center of 
the insulating layer of each MTJ. 



MgO. This could be due to the hybridization between 
O 2p z and Mg 2p Zl generating the energy gap along the 
A line. Furthermore, the lowest n of MgAl2 04 showed 
a fc x dependence similar to that for MgO. This suggests 
the possibility that the large TMR ratio stems from the 
coherent tunneling properties and the half-metallic char- 
acter of bcc Fe on the Ai state in Fe/MgAl 2 O 4 /Fe(001), 
like in Fe/MgO/Fe(001) MTJs. These results are incon- 
sistent with our calculation results, where the TMR ratio 
of the Fe/MgAl 2 O 4 /Fe(001) MTJ (160%) was one order 
of magnitude smaller than that of the Fe/MgO/Fe(001) 
MTJ (1600%). 

To further elucidate these results, we show in Fig. 5 
the band dispersion of bulk bcc Fe with a tetragonal unit 
cell (a=5.733A, c/a=0.5), corresponding to the unit cell 
of the electrode region in the Fe/MgAl 2 O 4 /Fe(001) MTJ. 
Since bcc Fe with a primitive lattice constant of 2.867 
Ahas minority-spin states with s orbital character at the 
boundary edge of the 2-D Brillouin zone of fen, these 
minority-spin states appear at fey =(0,0) because of the 
band folding that occurs when the in-plane cell size of bcc 
Fe is twice that of the primitive cell size. The new con- 
ductive channels at fey = (0,0) couple with the Ai evanes- 
cent states of MgAl 2 4 , giving the slowest decay in the 
barrier layer. This induces a relatively large conductance 
in the anti-parallel magnetization configuration, resulting 
in the reduction of the TMR ratio. Thus, we concluded 
that the TMR effect in Fe/MgAl 2 O 4 /Fe(001) MTJs is 
intrinsically smaller than that of Fe/MgO/Fe(001) MTJs 
because of the band folding effect in the minority-spin 
states of the Fe electrode. 

We calculated the projected density of states of 
the tunneling wavefunction (TDOS) on the Ai(s) 
and A 5 (p y ) states for the Fe/MgAl 2 O 4 /Fe(001) and 
Fe/MgO/Fe(001) MTJs. Figure 5 presents the majority- 
spin TDOS at the Fermi level at fey = (0,0) in the bar- 
rier regions as a function of the distance from the 
left junction of the MTJs in parallel magnetization. 
We found that the TDOS of the Fe/MgAl 2 O 4 /Fe(001) 
MTJ showed a slower decay for the wavefunction of 
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FIG. 5: (a) Majority- and (b) minority-spin band disper- 
sion of the bcc Fe electrode in Fe/MgAl 2 4 /Fe MTJs with 
a tetragonal unit cell (a=5.733A, c/a=0.5) along the [001] di- 
rection at /c// = (0,0), where the bold lines indicate the band 
with Ai symmetry in the tetragonal unit cell, which corre- 
sponds to the unit cell of MgAl20 4 . 



the Ai evanescent state than for the A5 state. The 
same was observed in the case of the Fe/MgO/Fe(001) 
MTJ. This is clear evidence of the significant contribu- 
tion of tunneling electrons with Ai symmetry, as com- 
pared to those with A5 symmetry in MgAl 2 4 . Fur- 
thermore, the decay rate of the TDOS of Ai states in 
the interior of the barrier layer is almost the same for 
both Fe/MgAl 2 O 4 /Fe(001) and Fe/MgO/Fe(001) MTJs. 
These results indicate that MTJs with a MgAl 2 4 bar- 
rier can potentially show a large TMR ratio, comparable 
to that for Fe/MgO/Fe(001) MTJs, if band folding of the 
Fe electrode can be suppressed. 

Finally, we discuss a possibility of the coherent tun- 
neling and the band folding in recent experiments of 
MTJs with the MgAl 2 4 barrier. Sukegawa, et al, 
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FIG. 6: (Color online) Majority-spin tunneling density of 
states projected on the Ai(s) and As(py) states at the Fermi 
level as a function of the distance from the left interface for a 
normal junction of (a) Fe/MgAl 2 O 4 (001) and Fe/MgO(001). 
The scale of the vertical axis is logarithmic. 



reported TMR ratios of 117% at RT and 165% at 
15 K for Fe/MgAl 2 O x (2nm)/Fe(001) MTJsI. They 
showed that the MgA^O^ barrier layer formed a 
normal-spinel structure with a (001) orientation, and 
the Fe/MgAl 2 O ;r (2nm)/Fe(001) MTJ had a epitax- 
ial relationship between the bcc-Fe(001)[110] and the 
MgAl 2 Ojc(OOI). This means that the coherent tunneling 
occurs through the evanescent Ai state in the barrier 
layer. Thus, the experimentally observed small TMR ra- 
tio in the Fe/MgAl 2 O ;r (2nm)/Fe(001) MTJ compared to 
that of the Fe/MgO/Fe(001) MTJi"^ can be attributed 
to appearance of new conductive channels in the anti- 
parallel magnetization due to the band folding effect. We 
believe that a modification of the periodic boundary con- 
dition along the in-plane direction by choosing ferromag- 
netic electrodes without the band folding at the junction 
with MgAl 2 04 is necessary to obtain the large TMR ratio 
in MTJs with the normal-spinel MgAl 2 04. 



MTJ, which was much smaller than that of the 
Fe/MgO(l nm)/Fe(001) MTJ (1600%). We concluded 
that the appearance of new conductive channels at 
fen =(0,0) by the band folding effect significantly reduced 
the TMR ratio of the Fe/MgAl 2 O 4 (lnm)/Fe(001) MTJ. 
The new folded channel in the minority-spin state of the 
Fe electrode coupled with the Ai evanescent state of the 
MgAl 2 04 barrier, which showed a slow decay in the bar- 
rier layer and contributed to the minority-spin conduc- 
tance. These results indicate that the TMR effect in 
Fe/MgAl 2 O 4 /Fe(001) MTJs is intrinsically smaller than 
that in Fe/MgO/Fe(001) MTJs. However, if the band 
folding effect can be suppressed by changing the peri- 
odic boundary condition along the in-plane direction, 
the Fe/MgAl 2 O 4 /Fe(001) MTJs show a large TMR ratio, 
comparable to that for Fe/MgO/Fe(001) MTJs. 
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